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Abstract: 

We study complexified Harmonic Oscillator with a position- dependent mass, termed as 
Complex Exotic Oscillator (CEO). The complexification induces a gauge invariance [T9~j 
ITT] . The role of PT-symmetry is discussed from the perspective of classical trajectories 
of CEO for real energy. Some trajectories of CEO are similar to those for the particle in 
a quartic potential in the complex domain [10J 
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Introduction: It came as a surprise when the early works [U El El 131 G2, El [7] showed 
that certain quantum theories with complex Hamiltonian have real spectra. Subse- 
quently this observation was explained 0, EH \TU\ from the fact that these systems en- 
joy the combined PT (parity and time reversal) symmetry. The consistency of these 
models as quantum systems was established [11] by constructing a positive definite in- 
ner product that generates unitary evolution. Later on there has been a lot of activity 
[T2l [T3l [HI [151 EE El EH] m the study of different aspects of PT-symmetric models. 

These models are referred to as "Crypto" -Hermitian models by Smilga [T9] . In [19] 
Smilga has also provided an alternative explanation to this behavior (of having real energy 
eigenvalues for a complex Hamiltonian) in terms of a gauge invariance. However, in an 
important earlier work by Mostafazadeh [2D] , it was observed in a general context that the 
real part of the Hamiltonian can generate the dynamics in a real phase space. In addition 
the imaginary part of the Hamiltonian, treated as a constraint, can generate symmetry 
transformations. The usage of certain class of coordinates in previous works [211 1221 E3] 
was also explained in [20] . 

In our analysis we shall follow the method developed in [T5] . The idea used in [19] 
is to complexify a real Hamiltonian system and subsequently treat the real part of the 
complex TC as the Hamiltonian H of the enlarged system. As a consequence the number 
of degrees of freedom of the new real Hamiltonian (H) is actually twice the original real 
Hamiltonian. By virtue of Cauchy-Riemann condition (for TC) and Hamiltonian equations 
of motion it is possible to show that both the real part H and the imaginary part G of TC 
(where TC = H + iG) , are separately conserved. Hence G acts as a First Class Constraint 
in the terminology of Hamiltonian constraint analysis of Dirac [23] and the presence of G 
ensures the equality of the degrees of freedom count of the Hamiltonian system before and 
after the complexification. In particular G = forces the energy to be real. This FCC 
is present in all such complexified systems and the gauge symmetry induced by it [23] is 
termed as Crypto-gauge symmetry [19] . In [19] it has been shown that specific features of 
some complexified models, (analyzed in terms of real variables), can be matched with their 



PT-symmetric counterpart in the complex plane [10J. The advantage of the formalism 
developed in [19] is that one can start with any real Hamiltonian model that is convenient 
and study its complex generalization. 

The study of Smilga [19] was generalized to more than one dimension by one of us 
in [25] for the complexified Harmonic Oscillator. It was shown that the straightforward 
generalization yielded a richer and qualitatively different constraint structure where both 
the First Class and Second Class Constraints (SCC) [21] are present. The presence of 
SCC induces a change in the symplectic structure. The additional constraints, besides 
G — 0, emerge from the demand that like energy, the angular momenta should also have 
real spectra. An interesting feature was revealed in this study [25]: the number of First 
Class and Second Class constraints are such that the degrees of freedom count before and 
after complexification remains unchanged. The present work with a different model - the 
Complex Exotic Oscillator (CEO) - (or equivalently complex Harmonic Oscillator with 
a position- dependent mass), also behaves in this way and it is clear that this feature is 
generic. Recall that Cauchy-Riemann condition (for 7i) and Hamiltonian equations of 
motion were all that were needed to show that both the real part H and the imaginary 
part G of 7i (where H = H + iG), were separately conserved and one can interpret G = 
as a constraint. In particular it was a First Class Constraint and induced the gauge 
invariance. However, we still have not been able to provide an analogous proof for the 
rest of the constraint algebra in higher dimensions. 

Exotic Oscillator (EO) - the parent model of the present study - has an interesting 
history. Similar models have been studied before in the guise of a Harmonic Oscillator 
with a position dependent mass [261 123 12EJ [291 ESS ED] . The model has also been quantized 
[271 [28j [29]. Our present work deals with the classical setup. We complexify the EO to 
obtain the CEO and show that in higher (three) dimensions the constraint structure is 
same as that of the Complex Harmonic Oscillator [25J . The major part of the paper deals 
with the CEO obtained from the one dimensional EO where we plot the trajectories of 
the path for fixed energy. It is found that the trajectories of the CEO are nontrivial 



generalization of the trajectories of the Complex Harmonic Oscillator [19] and one can 
interpret the latter as a special case of the former. Incidentally in some cases we find a close 
similarity between the trajectories presented here with those of the complex anharmonic 
(quartic) potential [32|l8ll9]fl 

We mention the role of PT-symmetry in our work and later provide a discussion on the 
possible classical analogue of Exceptional points [331 Ell ESI EE] in the present study. We 
follow the definition provided in [32j for PT-symmetry in classical mechanics in the context 
of our Hamiltonian system that yields a set of ordinary differential equations. Under PT- 
transformation, a generic complex function f(t) is transformed to /(—£)* such that for 
PT-symmetric functions f(t) = f(—t)*. Furthermore, an ordinary differential equation 
f(t) = F(f(t),f(t)) will be PT-symmetric if for every solution f(t) of this differential 
equation, f(—t)* is also a solution of the same equation. The same definition will hold 
for a system of ordinary differential equations. In general the solution of a differential 
equation may not respect the symmetry the differential equation enjoys. In case of PT- 
symmetric systems the solution will also be PT-symmetric if it remains invariant under 
reflection with respect to the imaginary axis. Now in the formalism used here [19] the 
variables we deal with are real and the complex conjugation of the solution does not play 
any role. However its role is taken up by the extra variables that appear in the process 
of complexification and PT-symmetry nature of a solution is judged in a similar way by 
noting if it is symmetric under reflection with respect to the axis related to the extra 
variable that comes from the complex sector. This will become clear as we proceed to 
draw the solutions or particle trajectories later in the paper. 

Although quantization of the classical model studied here remains outside the scope of 
our paper we would like to mention the following point: for quantization of a constrained 
system there is some amount of non-uniqueness involved since the two acts, i.e. imposition 
of the constraints and quantization are not commutative (see [19] for more details). 

The Exotic Oscillator and its Complexification: We start by describing the 
Exotic Oscillator (EO). The results presented here are not new and can be extracted 
3 We thank Professor Hook for pointing this to us. 



from the more general settings provided in [23, E7J EU [29] . We only reproduce it in the 
simplest model to highlight the contrasting behaviors of Harmonic and Exotic Oscillators. 
Consider the Hamiltonian, H = — b(xp) 2 ) where (xp) = XiPi, x 2 = XiX^ p 2 = 
PiVi'i i = 1) •■) n an d the particle mass is taken to be unity. It has the interesting dynamics 
Xi = — 26(|(p 2 — b(xp) 2 ))xi = —{2bH)xi where the Hamiltonian (H) itself appears in the 
place of the frequency parameter for a normal Harmonic Oscillator l. 

The above system is classically symmetric under the PT-transformations as defined in 
Introduction. The Hamiltonian can be obtained from the Lagrangian L = |i; 2 + 2 (i^bx 2 ) - 
A particular solution x« = AiSin(ut) obeying the condition bA 2 = 1 leads to the time 
independent Hamiltonian H = \uj 2 A 2 = 

Hence the above solution is bounded with the parameter space of Ai being a (hyper- 
) sphere A 2 = | but notice that the phase space motion of EO is unbounded since 
Pi = Co "( ut ) Ai . On the other hand a normal Harmonic Oscillator (HO) H = \{p 2 + ax 2 ) 
has a solution Xi = AiSin(u)t) for u 2 = a with H = \ijj 2 A 2 . Hence for HO there is no 
restriction on the amplitude A$ but the frequency u depends on the parameter 'a' present 
in the HO Hamiltonian. This is in contrast to the EO case (as shown above) where there 
is no restriction on frequency u but the amplitudes Ai are not arbitrary. It depends on 
the parameter '6' present in the EO Hamiltonian. Also, the phase space trajectories are 
bounded and unbounded for HO and EO respectively. 

In the present work we shall consider some aspects of phase space profile when we 
complexify the EO. This needs a stabilizing term for the EO and the simplest remedy is 
to introduce the Harmonic Oscillator potential with strength 'a', such that 

H= -{p 2 -b(xp) 2 + ax 2 ), (1) 

/Li 

with Xi = —[2bH + a(l — 2bx 2 )]xi. Now we can compare the relative effects of the two 
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For the sake of curiosity, we note that in general, for 



H = [\{p 2 - b{xpf)] v =h v ^Xi = -{2v 2 bh^-V) Xl 



for a constant 



interaction terms simply by tuning the parameters 'a' and l b\ Clearly this extended 
system also enjoys PT-symmetry. In the following analysis we deal with the Complex 
extension of EO and classify the constraints particularly for three dimensional case. 

We consider the CEO in the prescribed way [19], [25] by replacing x i: pi by Zi, 7Tj respec- 
tively, in the Hamiltonian ([!]), to yield 

H(n il z l ) = ^(n 2 -b(zn) 2 + az 2 ). (2) 

The identification of Zi,i\i is given in terms of real canonical degrees of freedom as = 
Xi + iyi ; 7Tj = pi — iqi and they satisfy the relations {xi,pj} = 5y ; {y^ qj} = The 
complex Hamiltonian 7i in (pQ) now reads, 

W(7Ti, z^ = H(p h qi, Xi, yi) + iG(p h qi, Xi, y { ) (3) 

where the real and imaginary parts are respectively, 

H = ^[(p 2 - q 2 ) + a(x 2 - y 2 ) - b[(xp) 2 + (yq) 2 - (xq) 2 - (yp) 2 + 2(xp)(yq) + 2(xq)(yp)] 



G = -(pq) + a(xy) - b[(xp)(yp) - (xp)(xq) + (yp)(yq) - (xq)(yq)]. (5) 

Customarily the real part H is taken as the Hamiltonian that generates the time evolution. 
The real and imaginary part satisfy the relation {G, H} = 0. We further restrict the 
system by imposing the constraint, 



G^O, 



(6) 



(The weak equality is interpreted in the sense of Dirac [21].) This simplifying choice is 
same as that of [19j but indeed one can (and should) consider more general models with 
G being some non-zero c-number. 



In a previous work by one of us [25] we considered the complex Harmonic oscillator 
and discussed the full constraint algebra by including the angular momentum as well. 
The angular momentum, given by Li = eijkXjPk, satisfy the relation = {Li,H} = 
in both the cases. Interestingly the rest of the algebra is also repeated here although the 
form of H and G are different. 

Lastly, it is worthwhile to study the equation of motion: 

Xi = -[a + 2bH - 2abx 2 + 2aby 2 ]x l + 2b[G - 2a{xy)]yi (7) 

Notice that for a = 0, that is without the Harmonic Oscillator potential term we obtain, 

Xi = -[2bH]x i + 2b[G]y i . 

This means that even for the CEO, the characteristic feature of the dynamics is preserved 
on the constraint surface where G — 0. 

Dynamics of one-dimensional Complex Exotic Oscillator (CEO): We study 
the CEO (JIJ) in one dimension, 

H=\ V \\-bx 2 ) + a -x\ (8) 

Obviously the system is considerably simplified. Also it is clear to see that it constitutes 
an oscillator with a position dependent mass. The equation of motion is 'x = —2[bH — 
a{bx 2 — \)}x. The PT-symmetry of the Hamiltonian is manifest. 

We now extend the system to complex domain with z = x + iy ;n — p — iq, 

n= l -^{i-bz 2 ) + a -z 2 . ( 9 ) 



The real and imaginary parts of 7i respectively are 

H = ^[a(x 2 - y 2 ) + (p 2 - q 2 ){l - b{x 2 - y 2 )) - Abxypq], 



(10) 



G = xy(a — b(p 2 — q 2 )) — pq(l — b(x 2 — y 2 )). 



(11) 



In one dimension G is the only constraint since {G, H} = and it is a First Class Con- 
straint inducing gauge invariance. We consider the case where G = 0. The Hamiltonian 
equations of motion are, 

x = p{l — b(x 2 — y 2 )) — 2bxyq, (12) 
p = bx(p 2 — q 2 ) + 2bypq — ax, (13) 
y = —q(l — b(x 2 — y 2 )) — 2bxyp, (14) 
q = —by(p 2 — q 2 ) + 2bxpq + ay. (15) 
For constant energy H = E the variables (x,p,y,q), must satisfy the conditions, 

E = ^[a(x 2 - y 2 ) + {p 2 - q 2 ){\ - b{x 2 - y 2 )) - Abxypq], (16) 

xy{a - b(p 2 - q 2 )) - pq{\ - b(x 2 - y 2 )) = 0. (17) 

These will help us in determining consistent initial conditions when we compute and 
sketch the trajectories (albeit numerically). 

Trajectories of one-dimensional Complex Exotic Oscillator: We now discuss 
some features of the trajectories. In the energy expression (fl6|) there are three parameters 
E, a and b. Another free parameter 'c' will appear from our choice of initial conditions. 
Throughout our analysis we will keep | E |= \ (since both positive and negative values 
of energy can be considered) and put a = 1 (that is strength of the stabilizing harmonic 
potential is fixed to unity). This will allow us to vary '6' (the strength of the exotic term 
(xp) 2 ) and 'c', the parameter that comes from fixing the initial conditions, mentioned 
above. Essentially 'c' signifies the freedom that we have in the complex domain so that 
we can choose initial conditions that allow the CEO to traverse paths which is forbidden 
for the EO, (the latter being restricted to the real space). This is the origin of the nested 



set of trajectories where in the core we have the real system confined to the real line in 
PT-symmetry framework [HI [9] (or real variable in the present case [IH]). To be more 
specific, changing 'c' is connected to a gauge transformation in the language of Smilga 
[19] that goes from one trajectory to another in the nested set of trajectories with the 
same physical parameters. We find that the trajectories are very sensitive to the initial 
conditions (i.e. the value of c). This feature has been stressed earlier in [37]. 

In Figure 1.1 and Figure 1.2, we study the phase space trajectories (x in abscissa and 
p in ordinate) of the real exotic oscillator (jSJ) , 

H = -p 2 (l - bx 2 ) + -x 2 = l = p 2 (l - bx 2 ) + x 2 
2 2 

for positive and negative values of b respectively. For our choice of parameters, '6' has an 
upper bound +1 for positive values. This is because for b = 1 the solution depicts isolated 
points p = ±1, x = ±1, and for larger positive '6' there are no closed orbits. On the other 
hand, there is no restriction on for negative values. For both Figurel.l and Figure 1.2, 
b = reproduce normal ellipse which we expect for the Harmonic Oscillator. In Figure 
1.1 the limit of '6' is < b < 1 and we find that the trajectories are diverging outside 
the normal ellipse, whereas in Figure 1.2, '6' goes up to b = —1500, and trajectories are 
converging inside the normal ellipse. Therefore the presence of the parameter '6' in the 
Hamiltonian is realized through the distortion in (x,p) phase diagram of the real exotic 
oscillator. 

Now we come to rest of the figures where the CEO is studied. From comparing the 
works of [H] (that deals with complexified models) and [H]] (that studies PT-symmetric 
models), it is clear that in our case the nature of the trajectories, as far as PT-symmetry 
is concerned, can be ascertained from the geometrical symmetry of the profiles. In all 
the figures we plot x in abscissa and y in ordinate following our convention z = x + iy. 
Hence, similar to [10] where real and complex parts of the coordinate were plotted in 
abscissa and ordinate respectively, trajectories that are PT-symmetric will be invariant 
under reflection about the ordinate. Note that this is same as the trajectories studied in 



PI as well. 

In Figure 2.1 we first reproduce the simple nested ellipses for the complexified Har- 
monic oscillator with b = 0, E = +| [13 [TOj . In Figure 2.2 similar example for CEO 
for positive non-zero b and E — +| are studied. It is clear that the non-zero exotic 
parameter '6' distorts the concentric ellipses. Trajectories for negative b are depicted in 
Figure 2.3. We notice a close similarity between Figure 2.3 and the anharmonic oscillator 
{H = \p 2 + x 4 ) in the complex plane studied in [32] . The graphs for negative energy is 
simply obtained by rotating the figures by |. In Figure 2.4 we fix the initial condition for 
E = +| to get one particular trajectory for given parameter values. The generic form of 
the initial condition (as referred by initial condition (A) in Figure 2) for these trajectories 
for positive and negative energies are 

x = c, p = 0, y — 0, q — 



x = c, p = 0, y = 0, q = 

respectively, where c is the free parameter we choose. All the trajectories with such initial 
condition are PT-symmetric closed orbits. It is interesting to note that the trajectory for 
b = —0.2499 in Figure 2.4 is identical to the limiting double cardioid depicted in [10] for 
H = p 2 + x 2 (ix) e , e = 2. 

In Figure 3 we choose a structurally different initial condition (as referred by initial 
condition (B) in Figure 3 and 4) of the form 

n, 1+46 
x = 1, p = c, y = 1, q = yc 2 + y^2> 

x=l,p = c, y = l, g=^_*Zi_ 

for positive and negative energies respectively. We notice that the orbits are not PT- 
symmetric and in general are not closed. However, for some specific choice of parameters 
(as in Figures 3.2) we do find closed orbits that are not FT-symmetric. This is inter- 
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esting because as mentioned in [32], (where examples of this type of orbits are shown 
for Harmonic Oscillator on complex domain), this occurrence is quite rare. On the other 
hand, Figure 3.3 is an example of an open orbit which is PT- symmetric in a restricted 
sense. Taken as a snapshot over several "rotations" one finds an overall left-right reflec- 
tion symmetry but since the orbit spirals inwards and outwards strict PT-symmetry is 
not maintained. Hence we feel that further studies are needed before one can make a 
general connection between closeness (openness) of trajectories with presence (absence) 
of PT-symmetry 

However, in the context of PT-symmetric quantum mechanics this type of broken 
{PT) symmetry situation is connected to the presence of Exceptional Points [331 El] 
where a singularity occurs in the parameter space due to the coalescence of two energy 
levels along with their wave functions. Note that this is distinct from a normal degeneracy 
where the wave functions of the degenerate levels are different. Characteristic features of 
the Exceptional Points survive in the classical trajectories of the corresponding classical 
problem [MJ [35[ [361 EH [38]. In fact Exceptional Points can be studied and observed 
in purely classical systems [39]. (Indeed it is all the more interesting if there exists a 
corresponding solvable quantum system.) In a classical model of two coupled damped 
oscillators this has been discussed by Heiss [34]. One can try to put our work in this 
perspective by interpreting our dynamical system ( fl2l to fT5l) as a very complicated non- 
linear extension of the two oscillator model of Heiss [34J. We provide an example of an 
Exceptional Point in Figures 4.1 and 4.2 where the values of b are different. Note that 
the nature and structure of the two trajectories differ drastically although the relevant 
points in parameter space are extremely close to each other. We interpret the point 
b = —0.18, c = 0.8 (in initial condition B) as an Exceptional Point since at this point the 
trajectories change from being strictly PT-symmetric (Figure 4.1) to non PT-symmetric 
(Figure 4.2). A difference between our system and the example studied by Heiss [31] 
is that the energy in our system does not become complex at any time since the choice 
of constraint G = is always present. As we have mentioned before, a more general 



choice for G may lead to interesting consequences. It would be worthwhile to pursue this 
problem further. 

Discussions: The discovery that one can replace the requirement of hermiticity to 
that of PT-symmetry in order to construct a Hamiltonian that supports unitary time 
evolution of quantum states and real energy eigenvalues has opened the possibility for 
considering various types of Hamiltonians that were rejected before on the grounds of non- 
hermiticity. Indeed, existence of PT-symmetry plays a role not only in models that are 
explicitly non-hermitian and violates either P or T symmetry individually, PT-symmetry 
is also relevant for models with Hamiltonians that are individually P or T symmetric. 
These models, such as Harmonic Oscillator or particle in a quartic potential exhibits 
closed trajectories that are PT-symmetric. It is rare that one finds [32] trajectories that 
are closed but not PT-symmetric. The complexified Exotic Oscillator that we have studied 
in the present article falls in this category. 

We have studied the complexified Harmonic Oscillator with a specific form of position- 
dependent mass. We have referred to it as Complex Exotic Oscillator. Our motivation 
for studying this particular form of effective-mass oscillator (in real space) is due to the 
fact that it is a PT-symmetric model (in the classical sense). It has an interesting form 
of Hamiltonian dynamics and furthermore is exactly solvable as a quantum system. 

We find that in general it is possible to have both open and closed trajectories for 
particles. The closed trajectories are mostly PT-symmetric but we have also found the 
example of closed orbits that are not PT-symmetric. We point out some possible connec- 
tion of this phenomenon with classical analogue of Exceptional points, that are associated 
with quantum mechanics when there is a crossover of eigenvalues from real to complex 
domain. Some of the PT- symmetric orbits resemble closely the orbits of a particle in a 
quartic potential extended in the complex domain [32] . 

Our aim is to study the quantized version [271 ES, EH] of the classical Exotic Oscillator 
in the context of Exceptional Point in order to see if the present findings can be corelated. 
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Collected Figure Captions: 

Figure 1: x vs. p plot: Distortion of the (x,p) phase diagram for positive and negative 
values of b and fixed E — \. In Fig (1.1) the limit of b is < b < 1 and in Fig (1.2) b is 
varying up to b = —1500. 

Figure 2: x vs. y plot: Family of trajectories for fixed E — \ and subjected to fixed 
initial condition (A); In Fig (2.1) value of b is zero; In Fig (2.2) b = 0.8; Fig (2.3) 
represents the family of trajectories for negative b (b = —10); Fig (2.4) is one particular 
trajectory for (A) and parameter values c = 2, b = —0.2499. 

Figure 3: x vs. y plot: Family of open (Fig 3.1) and closed (Fig 3.2) orbits lacking 
PT symmetry for fixed E — | and satisfying the initial condition (B). In Fig (3.1) 
c = 0.06, b = 0.5; In Fig (3.2) c = 0.06, b = 1; In Fig (3.3) b = 1, c = 1.2. 

Figure 4: x vs. y plot: Family of trajectories showing the Exceptional Point where 
the nature of the trajectory changes abruptly; In Fig (4.1) b = —0.09998, c = 0.8; In Fig 
(4.2) b = —0.18, c = 0.8. They both satisfy the fixed initial condition (B). 
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Figure 1: x vs. p plot: Distortion of the (x,p) phase diagram for positive and negative 
values of b and fixed E — -. In Fig (1.1) the limit of b is < b < 1 and in Fig (1.2) b is 
varying up to b = —1500. 
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Figure 2: x vs. y plot: Family of trajectories for fixed E — ^ and subjected to fixed initial 
condition (A); In Fig (2.1) value of b is zero; In Fig (2.2) b = 0.8; Fig (2.3) represents 
the family of trajectories for negative b (b = —10); Fig (2.4) is one particular trajectory 
for parameter values c = 2, b = —0.2499. 




Figure 3: x vs. y plot: Family of open (Fig 3.1) and closed (Fig 3.2) orbits lacking 
PT symmetry for fixed E = | and satisfying the initial condition (B). In Fig (3.1) 
c = 0.06, b = 0.5; In Fig (3.2) c = 0.06, b = 1; In Fig (3.3) b = 1, c = 1.2. 
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Figure 4: x vs. y plot: Family of trajectories showing the Exceptional Point where the 
nature of the trajectory changes abruptly; In Fig (4.1) b = —0.09998, c = 0.8; In Fig (4.2) 
b = —0.18, c = 0.8. They both satisfy the fixed initial condition (B). 



